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Abstract— This paper proposes a method for calculating the minimum amount of power load needed to shed and 
distributing it for each load bus in order to recover the frequency and voltage back to the allowable range. 
Based on the consideration of the primary control of the turbine governor and the reserve power of the 
generators for secondary control, the minimum amount of load shedding was calculated in order to recover the 
frequency of the power system. Computation and analysis of Voltage Sensitivity Index (VSI) of the load bus to 
prioritize distribution of the amount power load shedding at these positions. The lower the load bus have the 
Voltage Sensitivity Index (VSI), the higher the amount of load shedding will shed and vice versa. With this 
technique, frequency and voltage value are still within allowable range, and a large amount of load shedding 
could be avoided, hence, saved from economic losses, and customer service interruption. The effectiveness of the 
proposed method tested on the IEEE 37 bus 9 generators power system standard has demonstrated the 
effectiveness of this method. 

Keywords - Load shedding, Voltage Sensitivity Index (VSI), Frequency control, Primary control, Secondary 
control. 


The imbalance active power between the generation and 
the load demand causes a decrease the frequency in the 
power system. The monitoring and control system will 
immediately implement the control solutions to restore the 
frequency back to the allowable value, and improve the 
power system stability [1]. In [2], the primary and 
secondary control power plants are set by automatic 
controlled equipment or the power system operator. After 
implementing all possible control solutions that the 
system's frequency has not yet recovered to the allowable 
value, the most efficient method is to reduce the load [3]. 
The under-frequency load shedding relay (UFLS) is the 
traditional load shedding method used quite commonly in 
the current power system. In [4], the relays are set to 
operate whenever the frequency drops to a specified level 
and a fixed amount of load power is shed to restore the 
frequency. Using under frequency load shedding relay to 
disconnect the load bus will make in insufficient or 
excessive load shedding and take a long time to restore the 
frequency back to stable. This result will make damages 
for the suppliers and customers using the system's power. 
The authors in [5], [6] showed methods to estimate the 


I. INTRODUCTION 


amount load shedding, which based on the frequency 
reduction, or the rate of change of frequency (ROCOF). 
The combination of Intelligent load shedding methods has 
also been studied and developed such as Artificial Neural 
Network (ANN) [7], fuzzy logic algorithms [8], genetic 
algorithm (GA) [9] or particle swarm optimization (PSO) 
algorithm. These methods minimize the load shedding 
costs at steady state operation of pow system [10], [11]. 


A good load shedding program should be shed with the 
minimum number of load buses and as quickly as possible, 
and must meet the system's permissible frequency 
conditions. On the other hand, in large disturbances of the 
power system, frequency decay is often associated with 
voltage decay. Voltage decay at load buses reduces system 
load, therefore the reduction in frequency is slowed and 
the actual load shedding by UFLS is reduced relative to the 
level required [12] - [14]. In most practical conditions, the 
amount of load shedding may be more or less the amount 
of power the system needs to maintain the frequency, 
which may lead to cost losses as well as affects the objects 
that the system serves. 


This paper presents a new method of load shedding to 
minimize the amount of load shedding power. The load 
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shedding strategies based on the Voltage Sensitivity Index 
(VSI) to find out the priority and distribute the amount of 
load shedding power for each load bus. For load buses, the 
lower the load bus have the Voltage Sensitivity Index 
(VSI), the higher the amount of load shedding will shed 
and vice versa.er system [10], [11]. 


II. METHODOLOGY 

2.1 Overview the power systerm frequency responds 

The basic concepts of speed governing are best illustrated 
by considering an isolated generating unit supplying a 
local load as shown in Figure 1. 


Steam or 
Water 


Valve/Gate 


Pe 



Load Pl 


Fig. 1: Generator provides independent load 
The power system loads are a composite of a variety of 
electrical devices. For resistive loads, such as lighting and 
heating loads, the electrical power is independent of 
frequency. In the case of motor loads, such as fans and 
pumps, the electrical power changes with frequency due to 
changes in motor speed. The references [15], [16] showed 
the overall frequency-dependent characteristic of a 
composite load: 

A P e = A P L + DAco 

Nonfrequency-sensitive-load-change Frequency-sensitive-load-change 

where: A P L - is the load component does not depend on 


frequency, eg heat load, lighting, ...\DAw r - the change 

in load depends on the change of frequency, eg, motors, 
pumps, etcDeviation of power change;A^ r : Deviation 
of angle speed change; D: The percentage change in load 
with percentage of change in frequency varies, D is from 1 

-r 2%. 

The transfer function block diagram reflects the 
relationship between the load change and the frequency 
taking into account the governor characteristic, the prime 
mover and the load response shown in Figure 2 [17]. 



Fig. 2: The transfer function block diagram describes the 
relationship between the load changes and frequency 



Fig. 3: The effects of the governor speed droop and 
the frequency of load on the net frequency [18]. 

The transfer function relating the load 


change, A P L , to the frequency change, A CO , is 


Af(s) = AP L (s) 


-1 


Ms + D 


i + i 

R 


1 1 

1 1 

f 1 1 

K l + ST Gy 

y \ + sT CH y 

^Ms + D) 


( 2 ) 


Where: Kg the amplification stage; CO re j reference speed; 


Tch “charging time” time constant; A P Valve per unit change 

in valve position from nominal; Ms angular momentum of 
the machine in Laplace transform; R is equal to pu change 
in frequency divided by pu change in unit output; it is 
characteristic for the sliding speed 
adjustment; R = —A f / A P 

The purpose of system simulation in the form of a transfer 
function is to calculate the time response of the frequency 
deviation when the load change step is APl. From the 
above description, frequency deviation in steady state it 
means the values of the transfer function is determined for 
s = 0: 

The steady-state value of Af (v) may be found by: 


A/ steady state = lim[AA/(s)] = - 


-A PA — 


-A P T 


1 + 


( 3 ) 


R 


+ D 


When the power system has multiple generators with 
independent governors, the frequency deviation in steady 
state is calculated according to formula (4). 

-A P L (4) 


Af- 


11 1 r, 

-1-... H-h D 

R, R' R 
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Set 


P = 


— + 

R 


V 1 

D 


is the general frequency response 


K eq j 

characteristic of power system. It includes the adjustment 
characteristics of turbine mechanical power and load. 
From formula (4), obtain: 


Af = -AP L .J3 


(5) 


2.2 Primary and secondary frequency control in the 
power system 

Primary frequency control is an instantaneous adjustment 
process performed by a large number of generators with a 
turbine power control unit according to the frequency 
variation. Secondary frequency control is the subsequent 
adjustment of primary frequency control achieved through 
the AGC's effect (Automatic Generation Control) on a 
number of units specifically designed to restore the 
frequency back to its nominal value or otherwise, the 
frequency-adjusting effects are independent of the 
governor's response called the secondary frequency 
control. The process of the primary and secondary 
frequency control was shown in Figure 4. 



Fig. 4: The relationship between frequency deviation and 
output power deviation. 

Characteristic line (A) shows the effect of the governors: 
change the turbine power according to the change of 
frequency: 

In balance mode, the intersection of the generator 
characteristic line (A) with the frequency characteristic of 
the load line (D) determines the frequency fo equal 50Hz 
(or 60 Hz). When the load increases APl, the new 
characteristic line will be line (E): P t +AP. In addition, the 
intersection of the generator characteristic line (A) with the 
new load characteristic line (E) defines the new frequency 
fi. Here, fi< fo. Compared to the case where the generator 
does not have a governor, characteristic line (B), it is clear 
that: fi< fo. According to the characteristic line (A) of the 
generator unit, the governor does not prevent the 
frequency reduction: Af=fo-fi 


However, because the generator has the governor, it has 
helped to limit the large deviation of the frequency. 
Compared with the case the generators do not have a 
governor (B), the intersection with the new characteristic 
line of the load (E) will determine the frequency fi: 
f’i<fi<f 0 . Thus, the governor of the generator unit has the 
effect of reducing the large change of frequency known as 
the primary frequency controller. The efficiency of the 
primary frequency control depends on the slope of the 
speed-droop characteristic of the generator units. In the 
ideal case, the adjusting characteristic line (F) of the 
generator unit is vertical line; the frequency will not 
change until the power limit of the generator unit P n . 

The above characteristics of the primary adjustment 
process lead to the need for external intervention (by the 
automatic control device or by the power system 
operators) - that is the secondary frequency control 
process. The secondary adjustment characteristic is 
represented by the simultaneous shifting of the 
characteristic line (A) into the characteristic line (C) of the 
generator unit, with the slope unchanged. This adjustment 
is equivalent to the creation of a static vertical adjustment 
characteristic line (F). Thus, the secondary adjustment is 
within the rated power range of the generator unit to 
restore and maintain the frequency within the allowable 
value. 

2.3 Calculate the minimum load shedding power 
considering the control characteristic of turbine 
mechanical power and load 

In the 60Hz power system, the frequency deviation 
allowed Af p is 0.3 Hz (Af p = -0.3Hz). In the case of the 
power deficiency or imbalance between the source and the 
load causing the frequency difference, the frequency 
control will be implemented in the following order: 
primary control, secondary control. When the reserved 
power is used for secondary control and the frequency has 
not been restored to the permitted value, the load will be 
shed. Thus, from formula (5), the relationship between the 
permissible change in frequency, the amount of secondary 
control power and the minimum load shedding power 
PLSmin is calculated according to the proposed formula 
below: 


A f p P' [AP l (AFsecondary control 


( 6 ) 


In this Case, if (AP Secondary control APls min) ^ AP Sacondary max? 

then APLSmin=0, otherwise the minimum power load 
shedding is calculated by the formula below: 



-AR 


Secondary Max 


(7) 


Where: Af p is the permissible change in frequency (pu); 
PLSmin is the minimum amount of power required to shed 
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(pu); AP secondary control is the amount of secondary control 
power addition to the system. 

2.4 The Voltage Sensitivity Index (VSI) 

The main objective of calculating VSI is to find the most 
sensitive node of the system from voltage sensitivity point 
of view [19]. It considered a numerical solution, which 
helps operator to monitor how to shed the load to prevent 
frequency and voltage collapse. Nodes, having minimum 
voltage sensitivity index are selected and then, using 
equation (8) to calculate the voltage sensitivity index 
(VSI). 

VSI at bus i, is defined as [20]: 


VSL = 


Za-v,) ! 

1=1 


( 8 ) 


Where Vk is voltage at k th node and n is the number of 
nodes. 

After calculating the minimum amount of load shedding 
power required, the next issue determines the distribution 
the amount of load shedding power at load buses. The 
node with least VSI will be picked as the best location for 
the load shedding. That mean, the amount of load shedding 
at different load buses can be calculated in the same way 
as the principle of the load sharing in the parallel circuit. 
Due to parallel electrical circuits: 


Y 



1 1 1 


— + — + ... + — 
^1 ^2 


(9) 



p Z/U 2 z, 

Equivalent formula follows: 


( 10 ) 

( 11 ) 


P Z VST„ 

i _ ec l g e( l 

P ~ Z, VSL 


( 12 ) 


Here, the general formula calculates the amount of load 
shedding distribution at nodes according to the Voltage 


Sensitivity Index (VSI): P [s . 


VSIe, 

VSI , 


■Pr. 


(13) 


Where, n is the number of generator bus; i is the number of 
load bus; Plsi the amount of load shedding power for the i 
bus (MW); Pls min the minimum amount of load shedding 
power to the restore of frequency back to the allowable 
value (MW); VST: the Voltage Sensitivity Index (VSI) of 
the i bus; VSI eq : the equivalent Voltage Sensitivity Index 
(VSI) of all load buses. 

III. CASE STUDIES - SIMULATION AND RESULTS 


The effectiveness of the proposed method is tested on the 
IEEE 37 bus 9 generators system [21] which is shown in 
Figure 5. Total the active power and the reactive power of 
the system are 1024.31 MW and 215.94 MVAR 
respectively under normal operating conditions. The 
maximum active power and reactive power of the system 
are 1087 MW and 449 MVAR. The control solutions 
minimize the amount of load shedding and maintain 
steady-state frequency from 59.7 to 60 Hz. To test the 
effectiveness of the proposed method, the outage situations 
of the generator units are calculated, simulated and tested 
the parameters. In the case of calculations and simulations, 
the spinning reserved power to control the secondary 
frequency is also considered. All test cases are simulated 
on PowerWorld GSO 19 software. The results are 
compared with the results of the traditional load shedding 
method using under frequency load shedding relay. 

Apply the (5), (6), (7) formulas calculate the system 
frequency, the amount of primary and secondary control 
power and the amount of load to be shed. The results of 
the computation of the outage generator situations are 
shown in Table 1. 



Fig. 5: The IEEE 37 bus 9 generators test system. 


Table 1. The outage generators cases 


Name of 

Gen. Bus 

Freque 

ncy 

(Hz) 

In the 

allow 

range 

The 

primary 

control 

power 

value 

(MW) 

The 

second 

ary 

control 

power 

value 

(MW) 

The 

amount 

of load 

sheddi 

ng 

(MW) 

REDBUD 

59.97 

Yes 

10 

0 

0 

ELM345# 

1 

59.56 

No 

150 

12.5 

27.77 

ELM345# 

2 

59.56 

No 

150 

12.5 

27.77 

PEACH69 

59.62 

No 

112.31 

16.32 

11.89 

CEDAR69 

59.86 

Yes 

52 

0 

0 

BIRCH69 

59.79 

Yes 

80 

0 

0 
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PEAR138 

59.6 

No 

140 

17.38 

15.91 

PEAR69 

59.7 

Yes 

110 

0 

0 


In the test example, the sudden disconnection of the 
PEAR 138# 1 generator (bus 53) is simulated. Applying the 
equation (5) calculates the stable frequency value when the 
PEAR 138# 1 generator (bus 53) disconnects from the 
system. The frequency value is 59.57 Hz, and shows in 
Figure 6. 


60.1 



58.7 

0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 

Time (See) 

Fig. 6: The frequency of the system when the PEAR138#1 
generator disconnects 

After the PEAR138#1 generator suddenly disconnects, the 
frequency value is less than the allowable value. 
Therefore, the primary control and secondary frequency 
control which presented in section II.2 for frequency 
recovery should be implemented. 

The primary control process is done automatically by the 
turbine governor after the PEAR 13 8# 1 outage generator. 
The value of the primary control power of each generator 
turbine is shown in Table 2. 

Table 2. The value of the primary control power of the 


generators 


Generator 

The increased primary 
control power of each 
generator (MW) 

REDBUD69 (bus 14) 

5.2 

ELM345#1 (bus 28) 

23.6 

ELM345#2 (bus 28) 

23.6 

SLACK345 (bus 31) 

32.8 

PEACH69 (bus 44) 

21 

CEDAR69 (bus 48) 

7.5 

BIRCH69 (bus 50) 

11.2 

PEAR138 (bus 53) 

0 

PEAR69 (bus 54) 

15.1 


Total = 140 MW 


Because the recovery frequency is less than the allowable 
value, so the secondary control is implemented after the 
primary control. The spinning reversed power of the 
generators will be mobilized to perform the secondary 
control. In this case, the secondary control power is 17.38 


MW. The frequency of the system after the 
implementation of the secondary control is shown in 
Figure 6. 

Thus, after performing the secondary control, the recovery 
frequency is 59.65 Hz and has not yet returned to the 
allowed value. Therefore, the final solution is load 
shedding. Equation (7) is applied to calculate the minimum 
amount of load shedding power to recovery the frequency 
in allowable value. 


4/, 

P 


= 1.4- 


0.3 


0.00482x60 


-AP 


-0.1738 = 0.1591 pu = 15.91 MW 


After calculating the minimum load shedding power, the 
load shedding distribution at the load buses is calculated. 
The amount of load shedding at load buses based on the 
Voltage Sensitivity Index (VSI) value of all load buses. 
Calculation steps at section II.4 are applied to calculate the 
Voltage Sensitivity Index (VSI) value of all buses. The 
Voltage Sensitivity Index (VSI) value and the voltage at 
load bus (pu) at all load buses is shown in Figure 7. 

The priority load shedding distribution for each load bus is 
calculated based on the following principle: The lower 
VSI, the greater the amount of shedding power. Equation 
(7) in Section II.3 is applied to calculate the amount of 
disconnection power value at the load buses. 

1.04 



0.9 

0.02 


0.0195 



0.017 

But But But But But But But But But But Bu. But But But But But But But But But But But But But But But 
J 5 10 12 I) 14 19 16 17 18 19 20 21 24 27 JO iJ 94 37 44 48 90 99 94 99 96 

Name of Bus 

Fig.7. The VSI and the voltage at load buses (pu) after the 
outage generator (PEAR138) 

In order to compare the effectiveness of the proposed 
method, the load shedding method using under frequency 
load shedding relay is used to compare. The process of 
UFLS is implemented when the frequency reduces below 


www.ijaems.com 


Page | 622 






























International Journal of Advanced Engineering, Management and Science (IJAEMS) 
https ://dx. doi. org/10.22161/ijaems.512.2 


[Vol-5, Issue-12, Dec-2019] 
ISSN: 2454-1311 


the frequency setting threshold. The load is usually cut 
step-by-step based on the load shedding table that pre¬ 
designed based on the general rule and operator 
experience. These tables guide the amount of load that 
should be cut at each step depending on the decrease of 
frequency. These values are shown in Table 3. 


Table 3. The UFLS scheme using load shedding table[22] 


The 

steps 

UFLS 

Frequency 

(Hz) 

Time 

delay 

(s) 

The amount 

of load 

shedding (the 
percent of 

total load) 

(%) 

Total 

amount 

of load 
shedding 

(%> 

A 

59.7 

0.28 

9 

9 

B 

59.4 

0.28 

7 

16 

C 

59.1 

0.28 

7 

23 

D 

58.8 

0.28 

6 

29 

E 

58.5 

0.28 

5 

34 

F 

58.2 

0.28 

7 

41 

J 

59.4 

10 

5 

46 


The frequency and the rotor angle comparison between the 
proposed method and the UFLS method are presented in 
Figure 8. 


It can be seen that the proposed load shedding method has 
less the amount of shedding (76.28 MW) than the UFLS. 
Here, the recovery frequency value of the proposed 
method is lower than the UFLS method. However, this 
value is still within allowable parameter and acceptable 
range (59.7Hz). Especially, when considering the voltage 
value of the proposed method is equivalent to the UFLS 
method, although this method has less the amount of load 
shedding than UFLS method. This can be explained by the 
fact that a large load at load nodes with the lower VSI are 
disconnected causing the voltage to recover faster. Figure 
9 and Figure 10 show that the voltage at all nodes after 
load shedding according to the proposed method and the 
UFLS method is near the same. 


60.1 



Fig. 8: The frequency comparison between the proposed 
method and the traditional method 


The propose method 



0.5 

Bu» Bat Bat But But Bat Bat Bat But Bat But But Bat But Bat But Bat But But But But But But Bat But But But But Bat But Bat But Bat But But But But 

1 3 5 10 II 13 M 15 16 17 II 19 JO II J4 J7 I* 29 30 31 3: 33 34 35 37 31 39 40 41 44 47 4| SO S3 54 55 56 

Name of But 

Fig.9: The voltage at all bus when load shedding 
according to proposed method 

The UFLS method 
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Fig. 10: The voltage at all bus when load shedding 
according to the UFLS method 

IV. CONCLUSION 

A load shedding method considers to the primary and 
secondary control elements of the power plant to calculate 
the minimum amount of load shedding power and restore 
the frequency back to the allowable value. The proposed 
method ensures the frequency and voltage of the system in 
case of a severe generation-load mismatch. The selection 
of location and distribution of load shedding power at load 
buses are based on the Voltage Sensitivity Index (VSI) 
concept. The effectiveness of the proposed method has 
been demonstrated on a 9-machine, 37-bus system under 
test cases. The performance of this method is found to be 
better than that of a conventional UFLS scheme. The test 
results show that the proposed method results in reduced 
amount of load shedding while satisfying the operating 
conditions and limitations of the network. 
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